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Abstract 

We compute the running of the cosmological constant and Newton's constant taking into 
account the effect of quantum fields with any spin between and 2. We find that Newton's 
constant does not vary appreciably but the cosmological constant can change by many 
orders of magnitude when one goes from cosmological scales to typical elementary particle 
scales. In the extreme infrared, zero modes drive the cosmological constant to zero. 
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1. Introduction 

It is widely accepted that General Relativity (GR) can be regarded as a low energy effective theory and 
can be treated as a quantum field theory with an ultraviolet cutoff of the order of Planck's energy (see [1] 
and references therein). When energies of this magnitude are approached, corrections to Einstein's theory 
with terms quadratic or of higher order in the curvature will become relevant. Eventually, "new physics" 
is expected to set in. In the low energy limit the terms with the lowest number of derivatives of the fields 
dominate; assuming the effective action to be described by a local functional of the metric, the cosmological 
term will be the dominant one. Here one faces the problem of explaining why this term should be so small in 
natural units. For a general review of the problem see [2]. Several possible mechanisms involving quantum 
effects have been proposed [3,4,5]. One expects that the cosmological constant problem can be solved within 
the context of the effective theory, without reference to ultrahigh energy physics. 

In this paper we study the running of Newton's constant G = l/l&irn and of the cosmological constant 
A. We shall consider "matter" fields coupled minimally to a background metric; as usual in the background 
field method, the graviton (spin 2) field is treated like another matter field. All interactions between matter 
fields will be neglected, except for their interaction with the metric. The only self-interacting field will be 
the graviton, which in accordance with the discussion above is supposed to have an action consisting of a 
cosmological and Einstein term. 

We use the Wilsonian approach to the renormalization group, in a formulation given by Wetterich [6] and 
later applied in a gravitational context by Reuter [7] . In this approach one writes an exact renormalization 
group equation for a coarse-grained effective action, depending on some external energy scale k. This exact 
equation is then approximated by postulating a specific form for the action, containing only a few terms. 
In this way we obtain RG equations for n and A, containing contributions from all possible matter fields 
with spin up to two. These solutions can be solved in the high-energy (low curvature) and low energy 
(high curvature) approximations. In the high energy regime the behaviour of the couplings follows power 
laws dictated essentially by dimensional arguments. In the low energy regime the running of the couplings 
essentially stops, except for a residual effect of the zero modes, which drives the ratio A/k to zero. 

The plan of the paper is as follows: in section 2 we write down the exact equation, in section 3 we make an 
Ansatz for the effective action and derive RG equations for the cosmological constant and Newton's constant; 
in sections 4 and 5 we study the high- and low-energy limits respectively; in section 6 we extract gauge- 
independent information from our calculations. In section 7 we discuss the significance and the limitations 
of our results, and their relation to other works. 

2. The exact RG 

In a scalar field theory with quantum fields <f> propagating on a Euclidean background metric we 
define a modified generating functional for connected Green functions, Wk, by adding to the classical action 
a term ASk, depending on an external momentum scale k: 

e -w„U) = J W ) e - s W)-/i*-^(*) . (2.1) 

The new term is constructed so as to suppress the propagation of modes with momenta smaller than fc. We 
will refer to it as an infrared cutoff. It can be taken purely quadratic in the fields, with the general form: 

AS fc (0) = \j d 4 x^4>Rk (-V 2 ) <\> , (2.2) 

where Rk — k 2 F i s a suitable function of the Laplacian. For example one may consider exponential 

cutoff functions of the form 

where 

/*(2)=exp(-a(p) 6 ) . (2.4) 
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We allow b to have values between 1 (gaussian cutoff) and infinity (sharp cutoff). In general the results will 
depend on the form of the cutoff function, i.e. on the parameters a and b. 

In flat space, the effect of the addition of the term ASk is to replace in the propagators the momentum 
squared by the modified inverse propagator: 

q 2 + R k (q 2 )= 1 _ q p k{q2) =Pk(<l 2 ) ■ (2-5) 

The effective average action is defined as the Legendre transform of W k , minus the cutoff: 

T fe (0) = W k (j) - J d^xj^ - A5 fe (0) . (2.6) 

Hcuristically, one can interpret T k as a kind of coarse-grained action. The scale k can be understood 
physically as the resolution with which the system is observed by some apparatus. Modes with momenta 
larger than k cannot be directly observed and their effect is averaged out by the functional integral. The 
fact that certain results depend on the shape of the cutoff function just means that observing the system in 
different ways will show different pictures. Other results may turn out to be independent of the shape of the 
cutoff function; they are then said to be universal. 

One can easily prove that the average effective action satisfies the following exact renormalization group 
equation [6]: 

9 ^4 Tr (S + i?fe ) ^ 

This equation describes the flow of T k as the external parameter k is varied. It is important to note that 
due to the form of the cutoff term AS k , the trace is actually finite, so there is no need for any regularization 
or renormalization to make sense of it. 

The exact equation has much the same form as an improved one-loop renormalization group equation. 
In fact, at one loop the effective action would be 

r(i> = ilndet0 ; O = || (2.8) 

and therefore if we define a modified one-loop effective action by 

if ) = \ lndet(9 fc ; O k = = O + R k (2.9) 



we get 

d t T^ = ^(detOfcJ-^detOfc = ^Tr^ 1 ^ , (2.10) 



2' ' 2 

whence 



;>..r : = - i ( — _ + //,. I d t R k . (2.11) 



,(D 1„ ( 5 2 S 
k 2 



This is the same as (2.7), except for the replacement of the classical action S by at O -1 : this is the 
"renormalization group improvement" . 

The exact equations for fields with higher spin can be obtained following similar lines. The case of 
spin 1/2 fields was discussed in [8] for a flat background. In a curved background the Dirac operator is 
7 M V M , where V M V = d^tp + u^ b T a bip, with the spin connection given by c^ fc = e a x T^ lv e bv + e^c^e 6 * and 
Tab = | [7a, 7b]- The square of the Dirac operator is —V 2 + jR. For Dirac spinors in a gravitational 
background, writing AS k = J d A x ^R^^^W fj,)ip one obtains an exact RG equation of the form: 

( ' 5 2 V k {D) \ a j}{D) 



'"-^ ' "< ) « \ (2.12) 
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with the trace now running also on the Dirac indices. The same equation holds also for Weyl spinors, with 
the trace running over the appropriate subspace. In the case of Majorana spinors the RG equation has again 
the same form except for an overall factor 1/2 on the r.h.s. For massless spin 1 fields there is the extra 
complication of gauge invariance. This has been discussed in [9]. One can obtain an exact RG equation by 
interpreting S as the complete action, including gauge-fixing and ghost terms. For the photon field we will 
use a gauge fixing term of the form 

S GF (A,g) = ^J AVflfVA,) 2 , (2.13) 

with the corresponding ghost term 

S ghost (c,c,g) = J d A x V5c(-V 2 )c . (2.14) 

The resulting RG equation is 



1„ / 6 2 T 



a ' r '=in^ + N **-^ +fl H a,Rr '- (2 - i5) 

In the case of a massive (Proca) vector field one has to drop the ghost term. 

For nonchiral spin 3/2 (Rarita-Schwinger) fields, the RG equation takes the same form as for Dirac 
spinors except that there will be spacetime indices as well as Dirac indices. We write it formally: 

^ = - Tr (S + < S) )"^ 5) ' (2J6) 

Finally, the RG equation for a spin 2 field has been discussed in [7]. For diffcomorphism invariance we 
choose a background gauge fixing with background metric g pu . Denoting 7 the dynamical metric, 

Sgf(i,9) = 7^Z k J d 4 x^g^F p F v . (2.17) 

where F v = v 7 ^ (<5£V P — ^g pa V u ) (j pa — g pa ) and V denotes the covariant derivative with respect to the 
background metric. The constant Z k can be written as Z k = Kfc/ic, where n k is the scale-dependent Newton 
constant and R is the running Newton constant evaluated at some arbitrary reference energy scale k (see 
eqns.(3.3) and (3.19) below). The corresponding ghost term is 

S ghost (v, v,g) = J d A x Jg % (-g^W 2 - R^) v v (2.18) 

The exact equation has the form: 



-1 / ?9.t, \ -1 



(0± + R w\ a #» 



where 7^ = g plJ + Vhh^. 

Finally, we can consider a theory where all the fields mentioned so far are present at the same time. 
Assuming that the (quadratic) cutoff terms do not mix fields with different spins, the RG equation for the 
coupled system has the form 

d t T k — <9 t r fe | sca i ar + C?fT fc |spinor + 9 t T fc| V ector + 5*1^ |rs + d t T k Igravity (2.20) 

where the terms on the r.h.s. are given by equations (2.11), (2.12), (2.15), (2.16), (2.19). 
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3. The truncated gravitational RG 

The exact equation (2.20) describes a flow in the infinite dimensional space of all possible action functionals 
for the theory. It is therefore equivalent to infinitely many equations in infinitely many unknowns (the 
couplings of all possible terms in the Lagrangian). It is an impossible task to solve this equation. 

The simplest thing to do is to assume that the action has a certain form, parametrized by a finite number 
of coupling constants, plug it in the equation and solve the resulting finite system of PDE's. The virtue of 
this approach is that the truncated, finite system still contains genuine nonperturbative information. The 
flip side is that there is no control on the approximation. 

We separate the gravitational action, whose parameters we want to study, from the matter action, which 
has a fixed form and is not subject to RG flow: 

Tk = pmatter + pgravity _ ^ ^ 

For the matter sector we consider "free" matter fields, i.e. matter fields that are coupled minimally to the 
metric but otherwise have no other interactions. We will consider ns scalar fields cj>, no Dirac fields ip, 
riM Maxwell fields A^, np (Proca) massive vector fields and nps Rarita-Schwinger fields. Here the spinor 
fields are supposed massive and therefore non-chiral. The results for chiral spinors are obtained replacing 
no by n\y/2 {nw being the number of Weyl spinors) and urs by ucrs/^ (ucrs being the number of chiral 
Rarita-Schwinger fields). 

The action for these fields is taken to be: 

l -cj> (-V 2 + m 2 s ) <f> + i> ( 7 "V M + m D ) ^ + -F^F^ + ±. (V„A" f + c(-V 2 )c 

(3-2) 

the sum being over all particle species. The tensors F^ and G M „ are the field strengths for A^ and 
respectively and all covariant derivatives are calculated with respect to the metric 7. For the gravitational 
action we take 

rr Vity (7,5) = ^fc J d 4 x^(2A k -R( 1 )) + S GF (j,g) + S ghost (v,v,g) , (3.3) 

where n k — l/16irGh, A fc is the running cosmological constant, G k is the running Newton constant and the 
gauge fixing and ghost terms are as in (2.17) and (2.18). The background metric appears only in the gauge 
fixing term. 

As already remarked, the contributions of different fields are all decoupled. Given the explicit form 
of the action (3.2) and (3.3) we now have to evaluate the second derivatives that appear on the r.h.s. of 
equations (2.11), (2.12), (2.15), (2.16) and (2.19), specify the cutoff functions Rk and evaluate the traces. 

Since we want to isolate the coefficients of terms containing zero and one power of the curvature, it is 
enough to evaluate the r.h.s. for a background geometry of constant curvature, namely one satisfying: 

R^pa = —R(g^ P 9iya — g^gvp) ; R^v = -^Rg^ ■ (3-4) 

The scalar curvature R can be either positive or negative, corresponding to de Euclidean Sitter or anti-de 
Sitter space. We also assume that all matter fields have vanishing vacuum expectation value, and therefore 
vanishing energy-momentum tensor. With this Ansatz, the classical field equation coming from (3.1-2-3) 
just says that R = 4Afe. However, we shall stay off-shell as long as possible and therefore will not assume 
this equation for the background for the time being. 
For ns scalar fields we have simply: 

o p 1 n s d t Pk ,„ _x 

Od k Iscalar = "y lr (0) p + ^2 " (3.5) 



pmatter = £ j 
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The suffix (0) is a reminder that the trace is on a field of spin zero. 
For Dirac spinor fields we choose the cutoff term in the form 



AS k = J d 4 xdetei>R { k D) (^V^ 
where e M is the vierbein for the metric g and the function R k D ^ is 



(3.6) 



4 D Vv,) 




ft(-v 2 ) + f 



72 i E 



1 7 M V M . 



(3.7) 



This results in the replacement of the Dirac operator 7 M V M by -^===7^ 1 V AI , whose square is Pk + -j. (In 

flat space, q^ is replaced by yj Pk(q 2 )<ifj., where = g M /|g|.) Note that there are no ordering ambiguities in 
(3.7), as 7 M V M commutes with — V 2 on de Sitter space. * For no Dirac spinors we then have 



d t T 



n D 



t*- k | spinor 



Tr 



(1/2) 



d t P k 



P k + §+m 2 



(3. 



D 



the trace being over four-component Dirac spinors. For Weyl and Majorana spinors the remarks following 
eq.(2.12) apply. 

For fields of higher spin, carrying reducible representations of SO (A), it turns out to be convenient 
to decompose further each type of field into its irreducible components. Actually, on Euclidean de Sitter 
space, one can further decompose in representations of 50(5) [10,11]. For a vector field we use the 
decomposition = A^ + V M </>, where V^A^ = 0. In the functional integral this change of variables gives 
rise to a Jacobian determinant, which is J = ^/Dct (— V 2 ) [11]. For reasons of dimensional homogeneity it is 
convenient to define (f> = y/ (— V 2 )</>. The Jacobian for this further change of variables cancels the preceding 
one. t Taking into account also the gauge fixing term (2.13) we have for massless fields 



.(2) 



n M 
2 



d 4 x^/g 



^(-V 2 + f) 5 ^ + ^(-V 2 )^ 



whereas for massive fields 



^(2) 



n P r 
2 J 



d 4 Xy/g 



R 



m 



2 12 



(3-9) 



(3.10) 



In evaluating the trace that appears in Eq.(2.15) one has to pay special attention to zero modes. In fact, 
the operator —V 2 has a constant zero mode when acting on scalars which does not actually correspond to 
any eigenvalue of —V 2 on A^. This zero mode has therefore to be removed from the trace. Similarly, the 
ghost zero mode does not correspond to any gauge transformation of A^ and should also not be counted 
[4]. The cutoff is chosen as for scalar fields by requiring the replacement of —V 2 with Pk{— V 2 ). The ghost 
contribution partly cancels the longitudinal photon contribution, and we find 



dt^k I vector — 



■ 1 rn 9 t P k 

2 r(1) ^TT" 



+ nip 



(3.11) 



* Alternatively we could have chosen \J Pu (—V 2 + -R/4) in the numerator of (3.7). The square of the 
modified Dirac operator would then be Pk (—V 2 + i?/4). This operator commutes with 7^V M for any 
background field. 

t We observe that if we were to work with (f) instead, one should not introduce a cutoff term in the 
Jacobian determinant. This alternative procedure would lead to the same results as the method followed 
here. 
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where TV is a sum over all nonzero modes. 

For the Rarita-Schwinger field one has the decomposition: 



+ 77^ ; 7 m $m = 
*V =^+(v„-^V,)c; 



v*^ = o 



(3.12) 



with the corresponding functional Jacobian y Dot (—V 2 — j^R) [11]. Again it is convenient to define £ = 

V 2 — j^RC, so as to deal with fields all having the same (canonical) dimension and to cancel the functional 
Jacobian in the path integral. The linearized action is 



3 

+ 16 



C ( - ;; N , - 2m „.v ) C - 2C \/ -V 2 -^-^ (7 M V M - 2m KS ) V 



(3.13) 



The square of the operator acting on $^ is —V 2 + -f + "i|j g . The cutoff is chosen as follows: 
AS k -- 



^T R (RS)^T 



+ 



16 



Pk 



R 
12 



-v 2 - ^ U - 



(3.14) 



where 



i?f S) ( 7 ^V,) 



Pfe(-V 2 ) + f 



-V 2 + f 



1 7"V M . 



(3.15) 



Then vf> + AS k is equal to (3.13) except for the replacement of -V 2 by P fe (-V 2 ). The r.h.s. of (2.16) 
becomes 

nits m d t Pk 



dtFk\ns 



— Tr (3/2) p , , R 



P k + f + m RS 



(3.16) 



Finally, we come to the gravitons. As usual we define the graviton field ft M „ by 7 M „ = g M „ + 
The quantum fluctuation h^ u is decomposed into irreducible parts under SO(5): 

V = hl„ + V„£„ + V, y £ M + V„V v ct - i.9^V A W + J^ft , (3.17) 

where h — g^h^, £ M satisfies V a £a = and h£ v satisfies g^h^ = and V A ft^, = 0. 

We use instead of £ and a the variables £ M = ^ — V 2 — -j-R^ an d ^ = V— V 2 ^/— V 2 — |i?cr. As in the 
previous cases, the Jacobian for the functional change of variables from h^ u to h£ v , ct, ft, 



1/2 



J = Deti (-V 2 --i?J Dct ( -V 2 - -R) Det (-V 2 ) 



1/2 



2\V2 



(3.18) 



exactly cancels the Jacobian for the functional change of variables from ft^, , a, ft, to ft^, , a, ft. Thus 
there is no functional Jacobian to be taken into account when going from h^ v to h^ u , £ M , a, ft. 
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The quadratic approximation to the effective action in the neighbourhood of g, taking into account also 
the gauge fixing terms, has the form: 



,(2) 



— Ah 
2 



d'^gl t h T I _v 2 + ±R- 2A k ) h T ^ + I|„ 




a 



1 - 2a 



R - 2aA k ^ 



16 a 



& -V 2 - 



a 



a 



4a 



A, 



a 



k cr 



+ 2^—^dy/-^\ /-V 2 -^Rh- 
a — 3 



3a- 1 



h -V 2 - 



4a 



3a- 1' 



-A« 



3"" 3(3 -a) 

(3-19) 

Note that = k^/k plays the role of a wave function renormalization constant for the graviton. The 
running of Newton's constant can therefore be interpreted as due to an anomalous dimension of the graviton 
field. 

The ghost and antighost fields can be decomposed as the vector field A^; the ghost term (2.18) becomes 



.(2) _ 1 f 

k ~2j 



d 4 Xy/g 



R 



«J(-v 2 -y)rt + 0(-v 2 -^ 



(3.20) 



4" " n 2 
Note that we put the ghost wave function renormalization constant equal to one. 

We now have to decide the form of the cutoff terms ASk(g)- As before we assume that all the states 
with definite spin and parity, including the ghosts, have their propagators modified by the replacement of 
— V 2 with Pfe(— V 2 ). This is achieved by the following choice: 



1 



AS k = -Z k / d 4 x^\ (P k + V 2 ) h T ^ + -£„ (P k + V 2 ) £ 



a 



3 3 — a 
16 a 



1 



a - 



'Pk\ P 



l -R-^~^- 



V>--R h 



a (p k + v 2 ) a 

3a- 1 



3(3 - a) 



h (P k + V 2 ) h 



(3.21) 

so that Sk + AS'fc is identical to (3.19) except for the substitution of —V 2 by Pk everywhere. 

We now have to insert these formulae in the r.h.s. of (2.19). We decompose all the fields into normal 
modes of the operator —V 2 . Things go slightly differently depending on the sign of the curvature. For de 
Sitter space one has to take into account the fact that the first two normal modes of the field a do not 
contribute to h^ v and therefore have to be discarded (see also Appendix A). On the other hand even a 
constant h contributes to and therefore all the modes of h have to be included. In this way, the first two 
modes of h (i = 0, 1 in the terminology of Table A.I) remain singled out in the trace [12]. The ^-derivative 
acts on the explicit /c-dependence in P k and also on the running parameter Z k (we shall see a posteriori that 
this dependence is actually very weak and can be neglected). 

Putting all together and taking the traces leads to the following: 



d t T k \ 



gravity 



+ Vr 



d t P k 



(2) 



Pk 



^ Tr (D 



d t Pk 



jR-2A k 2 ^p je + is=lR-2aA k 



-Tr 



d t P k 



(0) 



d t Pk 




1 \ - 

+ 9 E 

Pk + V 2 



dtPk(\f } ) 



2 ^Pk-2K k ' 



Tr (1) 



d t Pk 



1. 



P^AT + 2 Tr '«; 



ft 

ft + v 2 



\R 



Tr', 



P k + - 2aA k 

d t P k 



**=±R - 2aA t 



E 



iti) P k -\R 

Pk(\f) - Ay 



(0) 
3 



0.1 Pki^f) 



4a 
3a-l 



At 



4a W (P k - 2A fc ) (P k + 



-R-2aA k ) 



(1 - 3a) ((3 - a)P k + (a - l)R/2) +4a(a + l)A fe )(P fe + V 2 ) 



3(1 - a?y/j\\[Pk - \R (VKyJPk -\P- v^V^-V 2 - ^Rj | 



(3.22) 



8 



where a double prime means that the modes i — 0, 1 arc omitted from the trace. In the case of anti-de Sitter 
space the spectrum is continuous; the traces will be over all normal modes and the terms involving isolated 
eigenmodes have to be dropped. 

The complete RG equation (2.20) is obtained by summing the r.h.s. of eqns (3.5,8,11,16,22). There is 
obviously no exact solution even to this truncated equation. In the next two sections we obtain approximate 
analytic results in the two limiting cases: k 2 ^> \R\ and k 2 <C \R\. Before doing that, however, we can 
say something on the general behaviour of the solutions, which can be guessed on the basis of dimensional 
considerations. 

We will assume from now on that we are not too far off shell, i.e. that is of the same order of 
magnitude as R. The most important contribution to the traces comes from the modes with frequencies of 
order k. This is because the functions d t Pk(z) and Pk(z) — z, decay exponentially when z ^> k 2 . In the 
region z w k 2 , Pk(z) w k 2 , so if we assume that either k 2 » \R\ or k 2 -C \R\, in the denominators in (3.22) 
we can neglect R and A& or Pk respectively. 

Similar considerations apply for the matter fields. In addition, in the matter contributions there are 
various possible regimes, depending on how the scale k compares to the masses of the matter fields. Generally 
speaking, if the mass of a matter field is greater than the scale k, that field decouples and gives a negligible 
contribution to the RG equation. This can be seen for example in the case of a scalar field: when k 2 -C m| 
we can neglect Pk in the denominator of (3.5). The remaining trace is proportional to k 2 for dimensional 
reasons, so the r.h.s. of (3.5) is of order fc 2 /m| <C 1. The situation is similar for the other fields. 

As k grows the RG equation passes through various mass thresholds, and every time this happens the 
contribution of a new field has to be taken into account. At a given energy scale k therefore one need 
only consider the contributions of fields with masses lower than k. Thus for very low energies one needs to 
consider only the massless fields will matter. 

4. The small curvature (high energy) limit 

In the case k 2 ^> R one is probing spacetime regions which are small relative to the curvature radius of 
the background metric. In such regions spacetime can be regarded as almost flat, with small corrections. 
One can therefore reliably use the heat kernel expansion to calculate the traces. As noted in the end of the 
preceding section, the dominant contribution to the traces comes from the region of momenta of order k, 
where Pk is itself of order k, so we shall expand the fractions in (3.5,8,11,16,22) to first order in R and A& 
and evaluate the coefficients with heat kernel methods. 

Let us begin by examining the contribution of the scalar fields. In the case ms <C k in the r.h.s. of 
(3.4) we can expand the denominator to first order in mi; there remains to evaluate Tr%^ and Tr^^. 

The trace of a function of —V 2 can be expanded in powers of the curvature, with coefficients given by the 
Mellin transform of the function itself (see Appendix A) . The result is 



d t T k : 



n s 1 



scalar 



2 (4 



Q2 



dtPk 
Pk 



dtPk 
Pk 



(4.1) 



where for any function / of the operator —V 2 , the integral Q n (f ) is defined by 



1 r°° 

Qn(f) = / dzZ n -\f(z) . 

T (n) Jo 



(4.2) 



The integrals appearing in (4.2) can be evaluated explicitly and are given in Appendix B. 

The contribution of the other fields can be computed in a similar manner. Whenever there are terms 
of order R in the denominator one has first to expand the fraction to first order in R, and then evaluate the 
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remaining traces using formula (A. 10). We have: 



d t r k \ 



(4^J d4xv ~ g { G 



^n s - 2n D +n M + ^n P - 4n RS ^j Q 2 (~p^~ 



+ (-7;n s m 2 s + 2n D m 2 D - ^n P m 2 P + 4n RS m 2 RS ^j Q 2 



d t P k 

p2 



+R 



1 1 1 1 \ n (d t P k 



,1 3 3 4 \ (d t P k 



(4.3) 



In a similar way one can calculate the gravitational contribution, which is 



d t T k 



gravity 



+ {6 + 4a)A k Q 2 —^r 



U RO ( 9tPk 
~ 24 RQ ' 



55 a. \un f dtPk 



+ 



d t Z k 



5Q 2 (^V(6 + 4 a) A fe 0/ Pfe + Vi 



1 



- t;RQi 



P k 
P k 



pi 



p2 



(4.4) 



where we have neglected terms of order RA k . 
Comparing (4.3) and (4.4) with the formula 



d t T k = J d*x^/g [2d t {n k A k ) ~ d t n k R] 



(4.5) 



we can read off the complete off shell beta functions for Newton's constant and for the cosmological constant. 
For Newton's constant 



d t K k = a x k 2 + r) k a 2 k 2 , 
where r] k = d t In Z k can be thought of as the anomalous dimension of the graviton and 



(4.6) 



1 



fli = 



a 2 



(4tt) 2 
1 



1 



1 



1 



1 



13 



12 US + 3 UD 24 " M " 8 HP + 24.' q{ ~ 



25 



1 



-n D + -n M + 



-n P - -n RS 



(4tt) 2 

Similarly for the cosmological constant 



55 
24 



a ] 1 2 



(4.7) 



d t A k = — ( cxk 4 + c 2 A k k 2 + ^d;m 2 fc 2 + %(c 3 fc 4 + c 4 A k k 2 ) 



(4- 
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where i = S, D, P, RS and 



Cl 

c 2 

C3 

c 4 
ds 



1 



(4tt) 2 
1 

1 5 
: (4tt) 2 2 
1 

: (4^)" 2 



1 1 3 1 , 

-ns - n D + -riM + ~ A np - 2n RS + - ) q 2 
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— rtn H riM 

3 24 



1 



13 



24 



1 



8 
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3«fiS + 7^ + a ) 92 



92 1 ; 



# + 1 



/47 
V 24 



(4^4^ 2 



(47T) 



2 n D<ll 



A 1 3 2 



(4^) : 



: 2n RS q 2 



(4.9) 

Let us find the simplest solution to these equations. Since we are assuming k 2 » \R\, with i? finite, we 
cannot use these off-shell equations for arbitrarily small k. So we will have to integrate them starting from 
some finite value k, where our variables have initial values R and A. 

In addition the effective theory we are using here is only valid much below Planck's mass, so we also 
assume R^> k 2 . We get as a solution 

k 2 ) . (4.10) 



1 



K + ^ai(/c 2 



,2 

Within this approximation % = < 1. Thus Kk does not change much as k runs over its allowed range. 
Since we are assuming that we are not too far off shell, |A^| w \R\ <C k 2 , so in (4.8) we can neglect the terms 
with coefficients c 2 and C4 with respect to those with coefficients c\ and C3. We then have 



4 R 



k 



(4.11) 



Remarkably, even though A fc has dimension of mass squared, it scales like k 4 . This is a consequence of the 
fact that the coefficient of R° in the Lagrangian is Hk^-k and Kk is almost constant. The consistency of the 
approximation requires that A -c k 2 . 

For a scalar field, the coefficient ci that determines the running of the cosmological constant is positive. 
This sign can be understood in terms of the Casimir effect: raising the IR cutoff removes normal modes 
and since the vacuum energy of bosonic oscillators is positive, it lowers the vacuum energy. For a rigourous 
argument in flat space in the context of the Wilsonian RG see [13]. The effect has opposite sign for fermions. 
The relative weight of the various fields in c\ is just their number of degrees of freedom (counted by 60) • 
Thus ci is proportional to the difference of the total numbers of bosonic and fermionic degrees of freedom. 
This may turn out to be positive or negative, depending on the theory. 

If Ci > 0, if we take k = and if we assume A — 0, we have an RG trajectory that lies entirely within 
the domain of validity of the approximation. On this trajectory A^ changes by a factor 10 240 as k varies 
between the inverse Planck time and the present value of the Hubble constant. Neighboring trajectories will 
have similar behaviour for sufficiently large k. We shall discuss these results further in section 7. 

5. The large curvature (low energy) limit 

In this section we assume R > 0. Because Euclidean de-Sitter space is a four dimensional sphere, in the 
limit k 2 <C R the scale probed is larger than the universe itself. The spectrum of the operators involved is 
discrete, and even the lowest eigenvalue is of order R (see Table A.I). When the IR cutoff is much smaller 
than R, varying it does not change anything, since all nonzero modes have already been taken into account. 
So the only possible effects must come from the zero modes. Since the functions f(z) appearing in the traces 
in (3.10) are exponentially decaying functions of z when z > k 2 , in the limit k 2 / R — > the traces are equal 
to /(0) if the operator has a zero mode, and zero if the operator does not have zero modes. The operator 
— V 2 on the sphere has zero eigenvalues only when acting on scalars, so the extreme low energy behaviour 
will be dominated by the scalar sector: a priori, we could have contributions from the scalar fields <fi, the 
longitudinal part of the photon, the electromagnetic ghost, the components a and h of the metric fluctuation 
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and the longitudinal part of the gravitational ghost. However, as noted in section 3, the zero mode should 
be removed from the trace over the longitudinal photon, the electromagnetic ghost and the longitudinal 
gravitational ghost. Thus we could only have contributions from genuine massless scalars and from the trace 
of the graviton. 

For z^Owe have: 

^)-S(f)'" 1 ( 1 + °(f)> <«W-£(FrH(£)') <"> 

When b > 1 these functions are divergent at z — 0. In this case the dominant terms in the scalar traces will 
be the ones containing the highest powers of P k . We thus find 



d t T k - b{n s + 1) + 0(rj) , (5.2) 

where hs is the number of massless scalars. * The solution Tk = T + b(hs + 1) In (|) has logarithmic infrared 
divergences, whose meaning we will discuss in the next section. 

The most difficult regime to study is the one where « Aj. If we start from k = and let k 2 grow, we 
have to use the expansion (A. 12) and take into account progressively higher modes. The first correction to 
(5.2) will come from the first nonzero modes of —V 2 . Using (3.5,8,11,16,22) and Table A.I and neglecting 77, 
it has the form: 



o t r k = 



35 d t P k {lR) 10 d t P k {\R) 5 d t P k {\R) 



2 P fe (§P) + \R- 2A k 2 P k (\R) + 2s=l R - 2aA k 2 P fe (Ip) - ^A 
d t P k (\R) r d t P k (\R) 40 d t P k {\R) 



_ Q 4 ■_ ^ 

■ ■ ■ I I — / I I I- — / * ill , , . , 

(5.3) 



P k {\R)-\R P k {\R)-\R 2 ™P k {\R) + \R + m 2 RS 



1 

+ 



1Q _ d t P k {\R) K d t P k (\R) 



P k (\R) + \R P k {±R) 



10 d t P k (\R) 

+ TT n P' 



2 "P fe (i#) + iii + m 2 , 



8 a t P fc (ii?) 5 d t P k {\R) 

T. n D T-pr-r. ; r-=— + 



P fc (ii?) + iP + m 2 D 2 °P fe (iP) + m| 

Each term in this sum is a function of two variables, A k /k 2 and R/k 2 . When P//c 2 becomes large each term 
tends exponentially to zero. Since the n-th eigenvalue is of order n 2 R, the contributions of higher modes are 
very effectively suppressed. We shall examine these contributions further in the next section. 

6. Going on shell 

One worrying aspect of our calculations is the dependence of the results on the gauge parameter a. One 
notes that even the leading coefficient a\, which appears in the beta function for k in the high energy regime, 
is a-dependent. Similarly, in the low energy regime, the contribution of the first nonzero mode (5.3) is 
a-dependent. This phenomenon is common to all calculations of off-shell effective actions, and is due to the 
fact that we do not compute physical quantities. Physical observables would turn out to be a-independent. 

There are several possible attitudes in this respect. One is to note that in practice the value of a\ 
hardly matters. If one compares the value of Newton's constant at a planetary scale (or for that matter at 
a cosmological scale) with its value at a scale of a few centimeters, as relevant for example in a Cavendish 
experiment, one finds AG/G w 10~ 66 ai, which is utterly negligible anyway. This observation however only 
avoids the question. A second, more satisfactory answer would be to reformulate the theory in a way that is 
gauge independent. For example one could compute the running of the Vilkoviski-de Witt effective action, 
which is gauge independent by construction. In de Sitter space, the Vilkoviski-dc Witt effective action is 
obtained by simply putting a = [4]. 



* The result is not continuous in the limit 6 — > 1. For a gaussian cutoff (b = 1, a = 1), Pfc(0) = & 2 /2 and 
dtP k (0) = k 2 are finite and we can neglect them with respect to R and A k . All terms in (3.22) are of order 
k 2 /R or k 2 /A k , and therefore are negligible. 
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A third attitude, that we will assume in what follows, is to try and compute the beta function of other, 
physically more meaningful quantities. To guess what the right quantities may be, we recall the well-known 
fact that in one loop calculations the on shell effective action is gauge independent. In the present context, 
going on shell means just putting R = 4A. Since the exact equation that we are using has the form of 
an improved one-loop equation, we anticipate that the same may happen here. If we consider the general 
truncated RG equation, the matter contribution is a-independent anyway and we see from (3.22) that in 
the leading terms of the gravitational contribution a cancels out on shell, while in the term involving 77 this 
does not generally happen. However in the high energy regime, if we put R = AK k in (4.4) we find that a 
cancels: 



^ 1 247r2 Jo^ fdtPk\ 13. _ (d t P k \ 19. (d t P k 

C*l fc|gravity,onshcll = ^ 2 < y~p^~ J g" A *:Vl \ ~p^~ J + -^ A kW2 y~pT~ 



2 V Pk J 6 1 Pi 



(6.1) 



In the low-energy regime, the zero mode term was a independent and in (5.3), the contribution of the 
first nonzero mode becomes a-independent on shell. 

What is the physical significance of these results? The RG equation can be thought of as a vectorfield 
in the space of all coupling constants of the theory, defining the change of these couplings under a variation 
of some external parameter (here k). Suppose we want to calculate the value of this vector at a point in 
the space of couplings. If we use the background field method, we should use a background that solves 
the (quantum-corrected) field equations at that point in the space of couplings. The field equations will 
generally depend on the coupling constants, and therefore one should use a different background at each 
point in parameter space. We call the resulting RG equation the on shell RG equation. In the present case 
the effective action is postulated to be of the form (3.2-3) and the equation of motion says that R = AK k 
(Note that n k does not appear in this equation. We can assume that n is actually a constant). If we now 
fix a value for A^ , the field equation determines R; this in turn gives the spectrum of the Laplacian and via 
the equations given before, the variation of A k . The solutions of the on shell RG equation are flow lines, 
along which the background field evolves together with the coupling constants. This on shell RG was applied 
before in a cosmological context [14]. 

The meaning of the on-shell RG equation can be further illuminated by the following observation. Using 
that the volume of the four-sphere with curvature scalar R is J d 4 x^fg = 24-^-, we get on shell 

5 t r fc | onshcll = J d 4 x^g-[2d t k k K k - 2d tKk A k ] = -48ir 2 d t (j^j (6.2) 

Therefore, the combination of the coupling constants whose beta function is a-independent is their dimen- 
sionless ratio F k = K k /A k . This fits with the observation that in pure gravity A and k by themselves are not 
measurable, but their ratio is [15]. 

Let us now review the solutions of the on-shell RG equation. If we compute the beta function of F k in 
the high energy regime we get from (4.6-9): 

-cifc 4 + (ai-C2)A fe /c 2 -^d i m l 2 fc 2 +r7 fe (-c 3 /c 4 + (a2-C4)A fe fc 2 )^ . (6.3) 

One observes again that a cancels out in the combinations ai — c 2 and a 2 — C4. This equation is still 
a-dependent, but only in the subleading order. Within the approximation of eqs.(4.10,ll) we have 

d t F k = -cJ^Fl (6.4) 

This gives back (4.11), which is therefore a trustworthy result. 
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In the low energy regime, we have from (5.2) and (6.2) 



(l + n s )b , , 

which yields 

Afe = l-^A ln i (6 ' 6) 

For k — > we then have A& — » 0. This running is very slow, however. 

Let us now examine the contribution of the first nonzero modes. Neglecting the contribution of massive 
particles, (5.3) becomes on shell: 



35 cWfAfc) 10 d t P k (A k ) 10 d t P k {A k ) 5 d t P k (±A k ) 
2P fe (§A fe ) + §A fc 2P fc (A fc )~A fe 2 nM P fe (A fc ) + A fc 2™ M P fe (§A fe ) 



8 d t P k (±A k ) 5. 9 t P fc (|A fc ) 
4 nH/ p fc (|A fe )+A fe + 2" 5 P fe (|A fe ) 



(6.7) 



Each fraction in this sum goes exponentially to zero for A k /k 2 > 1 and is equal to b for A k /k 2 = 0. The 
shape of the function is similar to a damped exponential for b = 1 and tends to a step function for large 
b. The i-th modes will give similar contributions, with the argument in P k and dtP k of order i 2 A k and the 
overall coefficient of order i 3 . One understands from here the qualitative behaviour of the solutions. Putting 
the cutoff at momentum k means summing the first i modes, with i 2 = k 2 / A k . Adding higher modes will 
pile up more of these functions, so for A k /k 2 < 1 the r.h.s. will grow fast with k, roughly like i 4 , while for 
A k /k 2 > 1 only the sero mode contribution remains, and it is independent of A k /k 2 . 

7. Discussion 

In this paper we have computed the scale dependence of the gravitational action as a function of an external 
parameter k, having the physical significance of resolution of the apparatus that is used to make the mea- 
surements. Features of the system with length scales smaller than fc _1 cannot be observed and are therefore 
integrated over. Their effect is summarized in the effective action T k , that describes the dynamics of the 
modes with wavelengths larger than fe _1 . 

In a gravitational context, the renormalization group is often seen as the behaviour of renormalizcd 
quantities under global (constant) rescalings of the metric [16]. We will now relate this point of view to the 
one used in this paper. Let us consider for the sake of generality an action functional S(g^ u , <f), A), where 
is a background metric, cf> are quantum fields and A are coupling constants. Under a global scaling of the 
metric g^ u — ► cj 2 g M „, the action is invariant provided all fields and coupling constants are scaled according 
to their canonical mass dimension 



S(w V, 0, A) = S(g^,w d *<f>, uj dx \) (7.1) 

In fact this relation can be taken as the definition of the canonical dimensions and d\. 

By following the definition of the average effective action given in the preceding section, one can easily 
prove the following equality: 

T k (u; 2 g^,^X) = T^g^^^^Lu^ X) (7.2) 

Note that k behaves like any coupling constant with dimension of mass. This formula, which follows simply 
from dimensional analysis, relates the behaviour of the effective action under changes of k to the behaviour 
under scalings of the metric. 

For u — > oo the distance between two given points in spacetime goes to infinity, so this limit is called 
the infrared limit and vice versa the limit u> — > is called the ultraviolet limit. On the other hand the 
Wilsonian renormalization group describes the behaviour of the theory as k is varied and one calls the limits 
k — > and k — > oo the infrared and ultraviolet limit respectively. One sees from eq. (6.2) that for example 
uo — > oo is equivalent (modulo rescalings of the fields and coupling constants) to the limit k — > oo and thus 
the terminology of what constitutes the ultraviolet and the infrared is reversed with respect to considerations 
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of scaling the metric. (This was implicitly stated in [14], eq.(23c) where a scalar held theory was shown to 
have an ultraviolet hxed point for £ = 1/6, while in the literature [16] this is called an infrared fixed point) 

We have already discussed the dependence of the results on the gauge parameter a. There is an 
additional dependence of the results on arbitrary parameters, namely on the parameters a and b that define 
the shape of the cutoff function (2.4). In this respect, one observes that the numerical value of all integrals 
(given in Appendix B) depends on the parameters a and b. The running of the couplings is therefore not 
universal. This is in accordance with the fact that Einstein gravity can be thought of as a spontancusly 
broken GL{A) gauge theory, with a symmetry breaking scale of the order of Planck's energy [17]. The 
rcnormalization group flow of coupling constants in spontaneously broken theories below the breaking scale 
is always non-universal. It is therefore reasonable to expect that gravity will exhibit universal properties 
only beyond the Planck scale. 

Even though we do not have a completely consistent theory of gravity in the ultraviolet, one can apply 
the methods of this paper to other effective theories of gravity that could be relevant in that regime. The 
running of Newton's constant was calculated before in a gauge theory of gravity with torsion and dilaton, 
containing terms quadratic in the curvature [18]. The results obtained here are in agreement with those of 
[18] in their domain of applicability. One thus gets a picture of a dilatonic theory in which the v.e.v. of the 
dilaton runs linearly above Planck's energy and settles to an almost constant value as the energy drops below 
the Planck threshold. This value can be identified with the (almost constant) Newton constant discussed in 
the present paper. We mention that the gravitational contributions to the running of a scalar self-coupling 
was computed using similar methods in [19]. 

Another intriguing possibility is that the effective action for gravity at low energies is non-local as 
suggested by several authors in the past [1,4,20]. In particular in [4] it was assumed that the effective action 
could contain a term of the form l^lnl/, where V is the four-dimensional volume. The effect of this term 
would be a highly nonlinear relation between the bare and effective cosmological constant. It is possible 
that our calculation, though restricted to local Lagrangians, may contain some hints of similar effects. More 
precisely, the running cosmological constant discussed here should be identified, for k — > 0, with A e gf. We 
also saw in Sect. 6 that in the IR limit terms like must appear. The existence of similar nonlocal terms 
could have profound consequences at the cosmological level. 
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Appendix A: Spectral geometry on the sphere 



The spectrum of the operator — V 2 on the sphere, acting on fields of definite spin, can be computed 
from the eigenvalues of the Casimir operator of the group SO(5) on irreducible representations [11]. Table 
I collects the results that are used in this paper: 
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TABLE A.I: Eigenvalues of —V 2 on the sphere 


Spin 


Eigenvalues 


Multiplicity 





\\ 0) = ±Ri(i + 3) ; *=0,1,... 


i(« + l)(« + 2)(2z + 3) 


1/2 


Af /2) = &R(i + I) 2 ; i=l,- 


f*(* + l)(* + 2) 


1 


Af } = ±R{i 2 + 3i - 1) ; i=l,... 


±i(z + 3)(2z + 3) 


3/2 


Aj 3/2) = ^i2(i + l) 2 ;i=2,... 


|(i-l)(i + l)(i + 3) 


2 


Af } = ±R(i 2 +3i-2) ; *=2,... 


f(*-l)(* + 4)(2* + 3) 



In this table the spin 1/2 and spin 3/2 fields are irreducible and therefore chiral. For Dirac spinors and 
nonchiral Rarita-Schwinger fields the multiplicities have to be doubled. Formally for a function of —V 2 wc 
have Tr/(— V 2 ) = J2i /(^») where Aj are the eigenvalues of —V 2 . The best studied such trace is the trace of 
the heat kernel e~ s (~ v ). For small s it has an expansion of the form: 

Tre -s(-V 2 ) = Bq( _ v 2 )s -2 + £ 2 (_ V 2 )s -l + B 4 (-W 2 ) + .... (Al) 

where B n = J d 4 x^/gtrb n and b n is a polynomial tensor involving the n/2-th power of the curvature. The 
coefficients of these polynomials depend on the space the operator is acting on. For operators acting on 
unrestricted scalar, vector, tensor and spinor fields, the coefficients can be computed using the formula [21]: 

(A2a) 

(A.2b) 

(A2c) 



bo = 


1 

(4tt) 2 1 


b 2 = 


1 R 


(4tt) 2 6 


b 4 = 


(4tt) 2 [( 



U80 - W*"- - + + lv,v^) 1 + 1 v 

where is defined as the commutator of covariant derivatives acting on the field in question. 

For fields subject to differential constraints, such as the spin 1, 3/2 and 2 fields Aj, and h^ v , 
the coefficients can be computed from the coefficients of the differentially unrestricted fields, using the 
decompositions (3.12), (3.17) and Table A.I. 

We give here a derivation appropriate to the Euclidean four-dimensional sphere. On a vector field A^, 
using the formula 

-V 2 V^ = V p (-V 2 - jR^j <j> (A3) 

one finds that 

Tre- s (- v2 V =Tre- s( - v ^\ A T +Tre- s( - v2 -iV\^ - e sR ^ (A4) 

The last term (which in the expansion (A.l) would appear as a —1 in the coefficient B 4 ) is due to the fact 
that the operator — V 2 — \R acting on scalars has a negative mode </>o=constant which does not correspond 
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to any normal mode of the operator —V 2 on vectors. From (A. 4) one can compute B n (— V 2 )|i. Similarly 
for the Rarita-Schwinger field using 

-V 2 (v„ - ^ V„) £ = (v M - (-V 2 - C (A5) 

one has 

Tre- s( - v2 )\*=Tre-<- v ^\*T+Tre- s (- v ^L+Tre- s (- v2 -iV\ c -8 (A6) 

The last term corresponds to the eight zero modes of —V 2 — on spinors that do not correspond to normal 
modes of —V 2 on $ M . For the graviton we use the formulae: 

-V 2 (V„&, + V„£„) = V M (-V 2 - + V „ (-V 2 - ^Rj ^ (A.7) 

and 

-V 2 ( V M V„ - ig^V 2 ) <r - fv M V„ - ig^V 2 ) f-V 2 - HjA a (A8) 



(A9) 



to obtain 

-rr e - s (- v2 )|^+Tr e - s (- v2 -A«)| 5 + T r e - s (- v2 )U+Tr e - s (- v2 -i«)| CT 
- e 2sR / 3 - he sR ' z - 10e sR / & 



The last term comes from the ten negative modes of the operator —V 2 — j^R on £, which are the Killing 
vectors of 50(5) and therefore do not correspond to normal modes of h^ v (V M £„ + V„£ M = 0); the second last 
term comes from the five negative modes of —V 2 (with eigenvalue — ^R) proportional to the functions x a , the 
coordinates of a flat space in which the sphere is embedded. These functions satisfy 4V Al V„£ a = c^V 2 ^; 
the third last term comes from the (constant) negative mode of —V 2 — |i? on scalars, (with eigenvalue 
— |i?). * In this way one can compute the relevant heat kernel coefficients for —V 2 acting on irreducible 
representations of 50(5): 



TABLE All: Heat kernel coefficients 



Spin 





1/2 


1 


3/2 


2 


tr&o 


1 


2 


3 


8 


5 


tr& 2 


\R 


\R 










For non-chiral (Dirac) spinors and Rarita-Schwinger fields, the coefficients in the third and fifth columns 
have to be doubled. 

For a function f(z), one finds the expansion 

Tr/(-V 2 ) = Bo(-V 2 )Q 2 (.f) + B 2 (-V 2 )Qi(/) + B 4 (-V 2 )Q (f) + ... (A10) 

* The presence of the constant terms in (A. 4), (A. 6) and (A. 9) can be confirmed by explicit calculations 
of the coefficient B4 of the l.h.s., and making a comparison to the zeta functions of the operators appearing 
on the r.h.s., whose spectra are known. 
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where 

Qn{f) = ==-r / dzZ n - X f{£) . (A. 

r <» Jo 

For large s the trace of the keat kernel is dominated by the lowest eigenvalues: 

Tre- s{ ~ v2) =e- sX " +er sXl + .... (A. 
and more generally for a function that decays fast enough for large s 

Tr/(-V 2 ) =/(Aq)+/(A x ) + .... (A. 

Appendix B: some integrals 

We collect here the result of some integrations which appear in section 4: 

ftflA ,.4 1,4 4 „ , 2 



dtPk\ , 2 2 7 2 2 ( 1 



Q 2 ^j = ^ = ^_r^| ; (B. 

*(^)-*«-*M* r G)' 

>fe + V 2 \ , 44 , 4 1 ^(2 



^ = ^—^j-T - ; (B. 



Pk J 2 (2a) 2 / b b \b 
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